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> ; Abstract 
O 

o 
o 

SU{2) lattice gauge theory with four flavors of quarks is studied at nonzero 

o , 

■ chemical potential fi and temperature T by computer simulation and Effective 

• Lagrangian techniques. Simulations are done on 8^, 8'^ x 4 and 12^ x 6 lat- 

• tices and the diquark condensate, chiral order parameter, Wilson line, fermion 



- 

X 



energy and number densities are measured. Simulations at a fixed, nonzero 



■ quark mass provide evidence for a tricritical point in the /i-T plane associated 

with diquark condensation. For low T, increasing /i takes the system through 
a line of second order phase transitions to a diquark condensed phase. Increas- 
ing T at high /i, the system passes through a line of first order transitions 
from the diquark phase to the quark-gluon plasma phase. Using Effective 
Lagrangians we estimate the position of the tricritical point and ascribe its 

existence to trilinear couplings that increase with // and T. 
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Recently there has been a resurgence of interest in QCD at nonzero chemical potential 
for quark number. Arguments based on instantons and phenomenological gap equations 
support the old idea that diquark condensation and a color superconductivity phase 
transition 0] will occur at a critical chemical potential of slightly less than one-third the 
proton's mass. Unfortunately, some of these arguments and calculations could prove to 
be misleading because a proper understanding of the diquark state alludes the field : since 
diquarks carry color in real QCD, a quantitative understanding of confinement and screening 
is needed to estimate the energies of the states which control the phases of the system. In 
addition, the brute force lattice simulation method does not yet have a reliable simulation 
algorithm for these environments in which the fermion determinant becomes complex. 

Since true QCD at finite quark-number chemical potential cannot be simulated with 
current methods, theorists have turned to simpler models which can be simulated. One of 
the more interesting is the color SU{2) version of QCD which addresses some of the issues 
of interest [Q], |Q. In this model diquarks do not carry color, so their condensation does 
not break color symmetry dynamically. The critical chemical potential is also expected to 
be one-half the mass of the lightest meson, the pion, because quarks and anti-quarks reside 
in equivalent representations of the SU{2) color group. Chiral Lagrangians can be used 
to study the diquark condensation transition in this model because the critical chemical 
potential vanishes in the chiral limit, and the model has a Goldstone realization of the 
spontaneously broken quark-number symmetry. Lattice simulations of the model are also 
possible because the fermion determinant is non-negative for all chemical potentials. One 
hopes that these developments will uncover generic phenomena that will also apply to QCD 
at nonzero chemical potential. 

Preliminary lattice simulations of the SU{2) model with four species of quarks have 
been presented at conferences recently 0, and it is the purpose of this article to present 
the continuation of that work with a focus on the system's phase diagram at nonzero n 
and T. Earlier work on this model at finite T and fi was performed by More detailed 
discussion of our zero temperature, finite /x simulations including the spectroscopy of the 
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light bosonic modes will be presented elsewhere 0. We choose to consider the system's 
phase diagram here because of its elementary, fundamental character and because these 
simulations produced an unanticipated result compared to what has been found previously 
in various models PJ^. For fixed quark mass, which insures us that the pion has a nonzero 
mass and chiral symmetry is explicitly broken, we find a line of transitions separating a phase 
with no diquark condensation from one with a diquark condensate. Along this line there is 
a tricritical point where the transition switches from being second order (and well described 
by mean field theory) at relatively low to a first order transition at an intermediate n 
value. At large fi the line of first order transitions separates the diquark condensate phase 
at low T from a quark-gluon phase at high T. The transition is clear in diquark observables 
as well as the Wilson Line. A schematic phase diagram is shown in Fig.l. (Note that 
we have not included the line of first order transitions starting from the /i = 0, finite 
temperature transition, which exists for m small enough). We shall see that this simulation 
result, the existence of a tricritical point, has a natural explanation in the context of chiral 



Lagrangians. Following the formalism of [ITU | we shall argue that trilinear couplings among 
the low lying boson fields of the Lagrangian become more significant as /i and T increase and 
they can cause the transition to become first order at a /i value in the vicinity of the results 
found in the simulation. By hindsight, this behaviour should not have come as a surprise. 
fi plays the role of a second 'temperature' in this theory in that it is a parameter that 
controls diquark condensation, but does not explicitly break the (quark-number) symmetry. 
The existence of two competing 'temperatures' is a characteristic of systems which exhibit 
tricritical behaviour. 



3 



Schematic Phase Diagram of Diquark Condensation in the Temperature-Chemical Potential Plane 
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FIG. 1. Schematic Phase Diagram in the T-^u Plane. The thin(thick) hne consists of sec- 
ond(first) order transitions. X labels the tricritical point. 



We begin with a discussion of the simulation method and the numerical results. 
The lattice action of the staggered fermion version of this theory is: 

-5/ = E {xiPif^) + m]x + ^A[x^r2X + XT2X^]} 

sites 



(1) 



where the chemical potential is introduced by multiplying links in the +t direction by 



and those in the —t direction by e^^ [|Tl[]. The diquark source term (Majorana mass term) 
is added to allow us to observe spontaneous breakdown of quark-number on a finite lattice. 
The parameter A and the usual mass term m control the amount of explicit symmetry 
breaking in the lattice action. We will be particularly interested in extrapolations to A ^ 
for a given m to produce an interesting, realistic physical situation. This is the case that 
has been studied analytically using effective Lagrangians at vanishing temperature T but 
nonvanishing chemical potential /i [^], 

Integrating out the fermion fields in Eq.l gives: 



pfaffian 



\t2 a 

-A^ \T2 



det(^U + A2) 



(2) 



where 
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A = Pifi) + m (3) 
Note that the pfaffian is strictly positive, so that we can use the hybrid molecular dy- 



namics 



T3[] method to simulate this theory using "noisy" fermions to take the square root, 
giving Nf = 4. 

For A = 0, m 7^ 0, /i 7^ we expect no spontaneous symmetry breaking for small /i. For n 
large enough (/x > according to most approaches including Chiral Perturbation Theory 

10| , |T2|| ) we expect spontaneous breakdown of quark number and one Goldstone boson - a 



scalar diquark. (The reader should consult for a full discussion of the symmetries of the 
lattice action, remarks about spectroscopy, Goldstone as well as pseudo-Goldstone bosons, 
and for early simulations of the 8 flavour theory at A = 0.) 

Now consider the simulation results for the Nf = 4 theory on 8^, 8^ x 4 and 12^ x 6 lattices, 
measuring the chiral and diquark condensates ((x"'"t2x)), the fermion number density, the 
Wilson/Polyakov line, etc. We are also engaged in larger scale simulations on 12^ x 24 
lattices, where, in addition, we are measuring all local scalar and pseudoscalar meson and 
diquark propagators (connected and disconnected). These simulations will be presented in 
our paper on zero temperature results 0]. Some of these results, especially for relatively 
large quark mass have been previewed elsewhere P]. 

First consider measurements on a 8^, 'zero temperature' lattice. We simulated the SU{2) 
model at a relatively strong coupling (3 = 1.0 to avoid finite size and temperature effects. 
(More extensive simulations at /5 = 1.5 and m = 0.1 also on an 8^ lattice will be presented 
in our zero temperature paper.) The quark mass was m = 0.05 and a series of simulations 
were done at A = 0.0025, 0.005, and 0.01 so that our results could be extrapolated to 
vanishing diquark source, A = 0. Simple linear and spline extrapolations were done and 
these procedures appeared to be sensible for these exploratory calculations. In the future 
more theoretical control over the A ^ limit will have to be developed. None of the 
conclusions to be drawn here will depend strongly on the limit. In fact, everything we say 
could be gathered from our data at a fixed (small) A value, 0.005 say. However, only in the 
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limit of vanishing A do we expect real diquark phase transitions (at least when the transitions 
are second order), so it is interesting and relevant to begin investigating the A ^ limit. 

Two further difficulties occur for small A. First, the Hybrid algorithm slows down as A 
is taken small and diquark observables are calculated. Analogous problems are familiar in 
QCD simulations at vanishing /i but small m because the Dirac operator of the standard 
lattice action becomes singular in the chiral limit |jl3[|. Even these exploratory simulations are 
quite CPU intensive for this reason, although they are well within the capabilities of standard 
PC's, and workstations. (Most of the results reported here were done on SVl workstations 
at NERSC and the T90 vector processor at NPACl.) Second, the Hybrid algorithm suffers 
from systematic errors proportional to the discrete time step used in integrating its stochastic 



differential equations forward in Monte Carlo time |jT3|. Most of the simulations reported 
here used a time step dt = 0.005 to control these errors. Future, accurate simulations on 
larger lattices will need even smaller time steps, dt = 0.005 proved adequate here but 
systematic errors are under investigation and will be reported on at a later date. 

In Fig. 2 we show the diquark condensate plotted against the coupling f3. The quark mass 
m was set at 0.05, the coupling was (3 = 1.0 and A has been extrapolated to zero using raw 
data at A = 0.0025 and 0.005. 

diquark condensate vs. mu , 8M, beta=1.0, m=0.05, lambda=0.000. 
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FIG. 2. Diquark Condensate vs. /x 
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We see good evidence to a quark-number violating second order phase transition in this 
figure. The dashed fine is a power law fit which picks out the critical chemical potential of 
fic = 0.2573(2). Since the coupling here is too large to reside in the gauge theory's scaling 
window and the lattice is too small to escape significant coarse-grain errors, we expect only 
semi-quanitative results. However, the power law fit is fair, its confidence level is 7 percent, 
and its critical index is Pmag = 0.32(5) which is somewhat smaller than the mean field 
result Pmag = 1/2, predicted by lowest order chiral perturbation theory [|l^, fl^. Note 
that the quark mass is fixed at m = 0.05 throughout this simulation, so chiral symmetry 
is explicitly broken and this transition is due to quark number breaking alone. In fact, the 
chiral order parameter (xx) is nonzero and varies smoothly from 1.00 to 0.50 over the critical 
region shown in Fig.2. The fermion number density, shown in Fig. 3, also shows the diquark 
continuous phase transition. The approximate linear dependence of the fermion number 
density with /i is the expected scaling behavior above the transition in lowest order chiral 
perturbation theory [|T^, p!2| . 

fermion number density vs. mu , 8M, beta=1.0, m=0.05, lambda=0.000. 
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FIG. 3. Fermion Number Density vs. /i 



Now we turn to the main focus of this letter, the phase diagram in the temperature- 
chemical potential plane. To begin, consider the small (8^ x 4) lattice results. 
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Consider a slice through the phase diagram at a fixed small temperature, for variable 
In Fig. 4 we show the diquark condensate for P — 1.3, m = 0.05, with the data extrapolated 
to A = as discussed above. As /j, increases we find that a second order phase transition to 
a diquark condensate appears at Hc — 0.2919(4). The dashed line fit in Fig.4 has the critical 
index P^nag = 0.50(15), in good agreement with mean field theory. The fit has a confidence 
level of 21 percent. The diquark transition is also clear in a plot of the fermion number 
density, as expected. However, the chiral condensate varies smoothly over the diquark 
critical region, with (xx) ranging from 0.80 to 0.50. The Wilson Line is also smooth, 
varying from 0.10 to 0.20, with no obvious signs of a transition. 

Diquark Condensate vs. Chemical Potential, 8^^3x4, beta=0.13. 
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FIG. 4. Diquark Condensate vs. ji. 

Simulations at j3 — 1.1 and 1.5 gave similar conclusions. At /3 — 1.1, there was a 
continuous diquark transition at = 0.2544(3) and &X, (5 — 1.5 there was a continuous 
diquark transition at = 0.2950(3). So, as the system is heated {j3 increases), a larger 
chemical potential is needed to order the system into a diquark condensate. 

Next, simulations were run at fixed // and variable (5. First, consider a relatively large 
and interesting fi — 0.40. We show the diquark condensate as a function of /5 for /x = 0.40 
in Fig. 5. There is a clear jump in the condensate at j3 = 1.55(5), suggesting a first order 
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phase transition. 



diquark condensate vs. beta , 8^^3x4, mu=0.40, m=0.05, lambda=0.000. 
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FIG. 5. Diquaxk Condensate vs. /? for /x = 0.40 

A discontinuous transition is also strongly suggested by the behavior of the Wilson Line 
shown in Fig. 6. 

Wilson line vs. beta , 8^^3x4, mu=0.40, m=0.05, lambda=0.000. 
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FIG. 6. Wilson Line vs. /3 for /x = 0.40 



If this transition at relatively large /x and T is really first order, then there should be 
a tricritical point separating the region of continuous transitions at smaller /i and T. We 
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strengthened the evidence for a first order transition at jj, — 0.40 by running the algorithm 

on a larger (12^ x 6) lattice. In Fig. 7, 8, and 9 we show the diquark condensate, the Wilson 

line and the chiral condensate on the 12^ x 6 lattice. 

diquark condensate vs. beta , 12^^3x6, mu=0.40, m=0.05, lambda=0.005. 




FIG. 7. Diquark Condensate vs. p. 



Wilson line vs. beta , 12^^3x6, mu=0.40, m=0.05, lambda=0.005. 




FIG. 8. Wilson Line vs. /3. 
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chiral condensate vs. beta , 12^^3x6, mu=0.40, m=0.05, lambda=0.005. 
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FIG. 9. Chiral Condensate vs. (5. 

In all figures we see evidence for a discontinuous transition at /5 = 1.97(2) separating 
a diquark condensate and a quark-gluon plasma phase. Note that these are figures of raw 
data are taken at A = 0.005. However, since the transition appears first order, it should 
extend to finite (although probably small) A. 

(3 = 1.97(2) lies in the scaling window of the SU{2) lattice gauge theory with four species 
of dynamical quarks, so we can compare our /Xc = 0.40 to the theory's spectroscopic mass 
scales ||14|- In fact, m.„/2 = 0.31(1), so finite T effects have apparently raised /ic somewhat 
from its zero T value, m^/2. Such effects will be discussed below in the context of Effective 
Lagrangians. 

The exact location of the tricritical point should be found in larger scale simulations. 
On our small 8^ x 4 lattices we have some evidence that it lies between fi = 0.30 and 0.40. 
In particular, we repeated the simulations that produced the data in Fig.'s 5 and 6, but at 
/i = 0.30 and produced Fig.'s 10 and 11. 
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diquark condensate vs. beta , 8^^3x4, mu=0.30, m=0.05, lambda=0.00O. 



0.45 

0.4 

0.35 

« 0.3 

I 0.25 
§ 

^ 0.2 
? 0.15 
0.1 
0.05 




0.9 





1 


1 1 


1 1 


- 


"■■■i.. 




- 


- 






- 




































1 


1 1 


'■• f T 



1 1.1 1.2 1.3 1.4 1.5 1.6 

beta 

FIG. 10. Diquark Condensate vs. /3 for /Lt = 0.30 
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FIG. 11. Wilson Line vs. /3 for = 0.30 



The dashed curve in Fig. 10 is a power law fit to the data which has a confidence level 
of 15 percent, a critical index Pj^ag — 0.59(15) and critical couphng Pc — 1-4740(5). Ac- 
cording to lowest order chiral perturbation theory, the critical behavior along the entire line 
below the tricritical point should be in one Universality Class, Mean Field Theory, while a 
tricritical point should have distinct indices, including = 1/4. It would be challenging 
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and particularly interesting to confirm the value /3^ag — 1/4 in the next round of lattice 
simulations, on larger lattices closer to the theory's continuum limit. 

Is this tricritical point suggested here expected theoretically? To answer this question 
we studied the possibility of a first order phase transition between the superfiuid and the 
normal phases using the Landau approach. Our goal was to analytically understand the 
emergence of the first order phase transition observed in the numerical simulations at high 
enough /i, not to describe the second order phase transition for chemical potentials just 
above /Zc = m.,,/2. We constructed the Landau free energy that describes the behavior of 
the diquark condensate for chemical potential greater than fic- Within this approach there 
is indeed a first order phase transition that separates the superfiuid and the normal phases 
for high enough fi. 

In order to construct the Landau free energy, we have first to determine the symmetries 



of the system [jT5[. The low-energy effective Lagrangian that has been used to describe 
QCD with two colors and quarks in the fundamental representation at zero temperature 
and nonzero chemical potential in is solely based on the symmetries of the QCD 

partition function. For small enough temperatures, the symmetries remain the same as at 
T = 0. Therefore the effective Lagrangian provides us with a natural basis to investigate 
the symmetries of the system in the superfiuid phase at nonzero temperature. We only 
need to use the substantial properties of the low-energy effective Lagrangian. The reader 



is invited to consult [[I^jITI] for a more detailed discussion of the effective Lagrangian itself. 



The effective Lagrangian at lowest order reads 

,rl 



^TrS^S^S^S + 2^Tr5Sl"9oS (4) 



-^^TrfSfi^S^ + BB) - m^Re TrME 



where F^^ is the pion decay constant, and m^r is the pion mass. The matrix B retains the 
symmetries of the chemical potential term in the QCD Lagrangian, and the matrix M retains 
the symmetries of the mass matrix. Finally the field S = UIIU'^ contains both the minimum 
S of (I) and the fluctuations around it, i.e. the Goldstone fields, U = exp^iiTaTa/ F^,) (Ta 

13 



are the generators of the Goldstone manifold). At T = when /i reaches /Zc = ''^7r/2, the 
minimum begins to rotate: S = cosaSc + sinaS^, where Sc corresponds to the quark- 
antiquark condensate, and corresponds to the diquark condensate. The angle of rotation 
is given by a = for /i < fic, and by cos a = m^/Aii^ for /i > /ic |10ip^ . For fi > the 
ground state of the system is a diquark condensation phase. The diquark condensate breaks 
the symmetry generated by the U{1) baryon charge. 

In the superfiuid phase {fi > fic) and at zero temperature, if we expand the effective 
Lagrangian @ in the fluctuations of S (i.e. in the Goldstone fields), we find that there are 
terms that are even in the Goldstone fields, and that there are terms that are odd in the 
Goldstone fields. These odd terms come from the linear derivative term in (^); they are 
absent in the normal phase. These odd terms are always present in the superfiuid phase, 
even at small temperature. They involve the long range excitations that carry the same 
quantum number as the diquark condensate. Hence the Landau free energy can be written 
as [0 



L = Atx' + Bx' + Cx\ (5) 

where x is the order parameter, and t = T/T^ — 1 is the reduced temperature. The or- 
der parameter of the superfiuid phase is the diquark condensate. The Landau free energy 
contains odd terms. They are not excluded by the symmetries of the system: as we saw 
above in the effective Lagrangian, these symmetries allow the presence of odd terms in the 
Goldstone fields that carry the same quantum numbers as the diquark condensate. This is 
very different from the situation at /i = where the symmetry excludes odd terms in the 
pion fields both for three colors and for two colors From standard studies of critical 



phenomena, these odd terms are known to be able to drive a first order phase transition. 
We have expanded L to O(x^) assuming that all the essential physics near the critical tem- 
perature Tc appears at this order. Since x = in the symmetric phase above the critical 
temperature, there cannot be any linear term in x, the coefficient of the quadratic term 
must be proportional to t, and A has to be positive. Furthermore C has to be positive so 
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that Ixl — ^ C)0 is not a minimum of the Landau free energy. 

Some properties of the coefficients A, B, and C in (^) can be deduced from the effective 
Lagrangian (^). Once we are in the superfluid phase, the pion mass is practically irrelevant. 
For chemical potentials not very far above fic, the orientation of the condensate is almost 
completely in the diquark direction. From the effective Lagrangian (^) expanded to fourth 
order in the Goldstone fields and neglecting the pion mass, we infer that the coefficient in 
the Landau free energy (H) are of the following form 

A = a/i^, a > 

B = bfi{t + 1) (6) 
C = c^?/Fl, oO. 

The temperature dependence of the coefficient of the cubic term B comes from the SqE 
term in the effective Lagrangian (^. This coefficient B has to vanish at zero temperature. 
Finally we get that the Landau free energy density is given by 

L = at^\^ + h^i{t + + c-^x'- (7) 

TT 

Since we neglected the pion mass, the chemical potential we are using here has to be 
understood as /x — yUc = — The Landau free energy (|^) can only be used to describe 

the transition from the superfluid phase to the normal phase. Therefore the Landau free 
energy (|^) is only valid for /i > /i^. For n — He ^ l^t = bFT^j^/ac, if t < {jJit/{l^ — l^c + 
^J{fx- fic)'^ - the minimum is at x 7^ 0, but for t > (/it/ (/i - /i^ + - /ij^ - /^?))^, 
the minimum jumps to x = 0. The Landau free energy (|^ therefore predicts a flrst order 
phase transition for /i > /ic + hF-^j^/ac. For fi — fic < /it the cubic term cannot drive a flrst 
order transition. From this analysis we thus expect that for fi — of the order of a 

flrst order phase transition separates the superfluid phase from the normal phase. 

At T = a second order phase transition separates the superfluid and the normal 
phases. From continuity arguments, the phase transition between these two phase has to 
remain second order for small enough temperatures. With the arguments given above, we 
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find that the phase transition is going to become first order above some chemical potential. 
We therefore find that there is a tricritical point in the phase diagram. Our analysis is 
too crude to tell the exact position of the tricritical point, but it predicts that it should 
correspond to fi — of the order of F.,^. 

This work will be followed up by several additional investigations. We certainly need to 
improve the numerical data by running on larger lattices, closer to the theory's continuum 
limit. Once that is done and more spectroscopy data is accumulated, we will be able to 
state our results in physical units and check many of the predictions of lowest order chiral 
perturbation theory quantitatively. It will also be interesting to study QCD with small 
chemical potentials associated with the three light quark flavors [|18[, ||19[. In particular, we 
expect the phase diagram at nonzero isospin chemical potential and temperature to be very 
similar to the one discussed here. Although we cannot attack the SU (3) theory with a large 
Baryon number chemical potential, these other situations can be studied both analytically 
and numerically, and interesting new phases of matter have been found there which should 
be investigated further. 

This work was partially supported by NSF under grant NSF-PHY96-05199 and by the 
U.S. Department of Energy under contract W-31-109-ENG-38. D.T. is supported in part 
by "Holderbank"-Stiftung. The simulations were done at NPACI and NERSC. B. Vander- 
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